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Holonomic Quantum Computation (HQC) is an all-geometrical approach to quantum infor- 
mation processing. In the HQC strategy information is encoded in degenerate eigen-spaces of a 
parametric family of Hamiltonians. The computational network of unitary quantum gates is real- 
ized by driving adiabatically the Hamiltonian parameters along loops in a control manifold. By 
properly designing such loops the non-trivial curvature of the underlying bundle geometry gives 
rise to unitary transformations i.e., holonomies that implement the desired unitary transformations. 
Conditions necessary for universal QC are stated in terms of the curvature associated to the non- 
abelian gauge potential (connection) over the control manifold. In view of their geometrical nature 
the holonomic gates are robust against several kind of perturbations and imperfections. This fact 
along with the adiabatic fashion in which gates are performed makes in principle HQC an appealing 
way towards universal fault-tolerant QC. 

PACS numbers: 03.67.Lx, 03.65.Bz 



> 
O 



r» : 
o . 
o ■ 
o ■ 

^ : 

Ok 

-i— > ■ 

G ■ 

s : 



I. INTRODUCTION 

It is by-now a generally accepted fact that the laws 
of quantum theory provide in principle a radically novel, 
and more powerful way to process information with re- 
spect to any conceivable classically operating device [Q. 
In the last few years a big deal of activity has been de- 
voted to devise and to implement schemes for taking ac- 
tual advantage from this extra quantum power. In partic- 
ular in Quantum Computation the states of a quantum 
system S are used for encoding information in such a 
way that the final state, obtained by the unitary time- 
evolution of S, encodes the solution of a given compu- 
tational problem. A system S with state-space TL (the 
Quantum Computer) supports universal QC if any uni- 
tary transformation U G U(7i) can be approximated with 
arbitrarily high accuracy by a sequence (the network) of 
simple unitaries (the gates) that the experimenter is sup- 
posed to be able to implement. The case in which S is a 
multi-partite system is the most important one as it al- 
lows for entanglement, a unique quantum feature that is 
generally believed to be one of the crucial elements from 
which polynomial or exponential speed-up occurs H. 

In the above picture of QC the realization of the quan- 
tum network is achieved at the physical level by turning 
on and off external fields coupled to S as well as local 
interactions among the subsystems of S. In other words 
the experimenter "owns" a basic set of time-dependent 
Hamiltonians that she/he activates at will to perform 
the suitable sequences of quantum logic gates. At vari- 
ance with such a standard dynamical view of QC more 
recently several authors considered geometrical and topo- 



logical approaches [|3|-|6j . The peculiarity of these propos- 
als is somehow striking: over the manifold C of quan- 
tum codewords one can have a trivial Hamiltonian e.g., 
H\c = 0, nevertheless, obtaining a non-trivial quantum 
evolution due to the existence of an underlying geomet- 
rical/topological global structure. The quantum gates - 
or a part of them - are then realized in terms of op- 
erations having a purely geometrical/topological nature. 
Besides being conceptually intriguing on their own, these 
schemes have some built-in fault-tolerant features. This 
latter attractive characteristic stems from the fact that 
some topological as well as geometrical quantities are in- 
herently stable against local perturbations. This in turn 
allows for Quantum Information processing inherently 
stable against special classes of computational errors. 

In this paper we shall give a detailed account of Holo- 
nomic Quantum Computation (HQC) introduced in Ref. 
§ and further developed in Refs. f§, §, (l(J and jS). 
In this novel gauge-theoretic framework one is supposed 
to be able to control a set of parameters A G M., on 
which depends an iso-degenerate family T of quantum 
Hamiltonians {H(X)}. Information is encoded in an n- 
dimensional eigen-space C of a distinguished H(Xq) G T . 
Universal QC |l3| over C can be then obtained by adi- 
abatically driving the control parameters along suitable 
loops 7 rooted at Ao- The key physical ingredient is pro- 
vided by the appearance in such quantum evolutions of 
non-Abelian geometric contributions jl4] U 1 G U (n) (n > 
1) given by holonomies associated with a w(n)-valued 
gauge potential A |r^ , |l5| . In other words quantum com- 
putation in the HQC approach is nothing but the paral- 
lel transport of states in C realized by the connection A. 
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Therefore the computational power in the HQC approach 
relies on the non-triviality of the geometry of the bundle 
of eigen-spaces of T over the manifold of control parame- 
ters, M. It must be noticed that also in Refs. even 
though the evolutions are not adiabatic, holonomies play 
an important role. Nevertheless the connections involved 
in those approaches are Abelian giving rise to [/(l)-valued 
holonomies e.g., Berry phases fli6fl . It is then clear that in 
order to achieve universality these geometric gates must 
be supplemented by standard dynamical operations. On 
the other hand in the HQC approach the whole quan- 
tum network is built by means of holonomies. In this 
sense HQC is fully geometrical. It is worth observing 
that the computational subspace C can be thought of as 
the lowest-energy manifold of a highly symmetric quan- 
tum system; from this point of view HQC is a kind of 
ground-state computation. This last remark points out 
the potential existence of a fault-tolerant |l7j feature of 
HQC due to energy gaps and even spontaneous relax- 
ation mechanisms. Further fault-tolerant characteristics, 
of HQC models we considered, are related to the fact 
that the holonomies J7 7 realizing quantum computations 
turn out to depend just on the areas of the surfaces that 
the generating loops 7 span on certain two-dimensional 
sub-manifolds. When this area is given one can consider 
even very large i.e., "far" from the identity deformations 
of 7, but as long as they are area-preserving no errors 
are induced. Moreover as far as the adiabaticity condi- 
tion holds U 1 does not depend on the rate at which the 
control loops are driven. Hence, even with respect the 
issue of timing, HQC is robust. 

The paper is organized as follows. In Section II the 
general theory of holonomic quantum evolution is pre- 
sented as well as its application to quantum information 
processing. A new pedagogical proof of the adiabatic the- 
orem is given and the properties of the holonomies are 
analyzed, facilitating their application to quantum com- 
puting. In Section III the CP™ model is presented as 
a theoretical holonomic arena allowing for the analytic 
evaluation of the connection A, its field strength F and 
a complete set of calculated holonomies. The mathemat- 
ical steps which enable such a calculation are given in 
detail providing an analytical method for calculating the 
holonomies (Wilson loops) much used in many areas of 
theoretical physics. In Section IV a physical model is 
presented based on quantum optics. Known optical de- 
vices as displacers, squeezers and interferometers are em- 
ployed as control devices performing coherent evolution 
of the laser photon states. The laser beams propagate 
in a Kerr medium which provides the desired degener- 
acy. Further similar applications to harmonic oscillator 
setups are discussed, which are experimentally more vi- 
able. In the Appendix a more mathematical approach to 
the Holonomic evolutions is presented. 



II. GENERAL THEORY 

In this section the general theoretical framework of 
HQC is reviewed. While the exposition relies partly on 
Refs. some proofs have been added which clarify 

the physical concept of holonomic evolutions making this 
subject more approachable to the quantum information 
community. 

A. Quantum Evolutions 

Let us suppose that we have at disposal a family T of 
Hamiltonians that we can turn on and off in order to let 
an TV-dimensional quantum system to evolve in a control- 
lable way. Formally, we assume T :— {H(X)}\ e M to be a 
continuous family of Hcrmitian operators over the state- 
space 7i = C^. The parameters A on which the elements 
of T depend will be referred to as control parameters and 
their manifold A4 as the control manifold, thought to be 
embedded in \R N . Indeed, one has 

N 2 

H(X) = i^2^ a (X)T a eu(N) 
0=1 

where the T a 's constitute a basis of the iV 2 -dimensional 
Lie-algebra u(N) of anti-hermitian matrices, and 
<i>: TVf 1 — i- u(N) is a smooth mapping that associates to 
any A in the control manifold a vector in u(N) with T- 
components ($i(A), . . . , <&at2(A)). 

The evolution of the quantum system is thought of 
as actively driven by the parameters A, over which the 
experimenter is assumed to have direct access and con- 
trollability. Suppose we are able to drive by a dynami- 
cal control process the parameter configuration A G M. 
through a path 7: [0, T] — > A4. Hence, a one-parameter 
i.e., time-dependent family 

:= {Hit) := [$ o 7 (t)]:i e [0, T]} c T (1) 

is defined. Notice that even the converse is true: any 
smooth family {-ff (^)}te[o.T] defines a path in A4 = 

ft . The quantum evolution associated to the time- 
dependent family (|l|) is described by the time-dependent 
Schrodinger equation id t \ip(t)) = H(t) \ip(t)) and hence 
it has the operator form 

U-y := T exp{-i f dtH(t)}eU(N) (2) 
Jo 

where T denotes chronological ordering. The time- 
dependent quantum evolution (j^), for a given map $, 
depends in general on the path 7 and not just on the 
curve 7 ([0, T]) i.e., the image of 7 in the control mani- 
fold. In other words the unitary transformation (^|) con- 
tains a dynamical as well as a geometrical contribution, 
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the former depends even on the rate at which 7 ([0, T]) is 
traveled along whereas the latter depends merely on the 
geometrical characteristics of the curve. 

From the physical point of view the parameters A rep- 
resent in general external fields and, for multi-partite 
systems, couplings among the various subsystems. To 
illustrate this point let us consider H := (C 2 )® N ^ C 2 ™ 
i.e., a iV-qubit system. Then a basis for u(2 ) is pro- 
vided by the tensor products T a :— ^fL 1 a ct( where 
a: {1, . . . , N} h-> {0, 1, 2, 3} and <t := 1, a\ := a x , & 2 '■= 
(T y , (73 := a z are the Pauli matrices. It is then clear that 
any a which takes a non-zero value more than once e.g., 
on otj 7^ describes a non-trivial interaction which gener- 
ates entanglement between the qubits i and j. Therefore 
the ability to manipulate the weight of the contribution 
of TVs in the decomposition of H(X), amounts to the ca- 
pacity of dynamically controlling many-body couplings. 
This goal is, of course, even conceptually more difficult 
to achieve than the control of the real external fields, 
namely the interaction associated to single subsystem 
generators T a . Finally, we stress that there is still another 
possibility; the control parameters A could represent on 
their own quantum-degrees of freedom e.g., nuclear co- 
ordinates in the adiabatic approximation for molecular 
systems, treated in some quasi-classical fashion. This sit- 
uation arises when one performs an adiabatic decoupling 
between "fast" and "slow" degrees of freedom, getting 
for the former a Hamiltonian that depends parametri- 
cally on the latter |l6| . In this case the control manifold 
M. is nothing but the classical configuration manifold as- 
sociated with a quantum system. 

Within this framework the requirements for imple- 
menting universal QC Jl3[ can be expressed in terms of 
the availability of paths. Universality is the experimen- 
tal capability of driving the control parameters along a 
minimal set { T j } f= 1 °f paths which generate the basic 
unitary transformations J7 74 's, i.e. the gates. By suf- 
ficiency of this set we mean the ability to approximate 
any U € U(N) with arbitrarily high accuracy by means 
of path sequences. 



B. Holonomies 

Now we recall some basic facts about quantum 
holonomies. A more mathematical approach can be 
found in Appendix A, where some by-now standard ma- 
terial has been collected aiming to make the paper as 
much as possible self-contained. 

The non-Abelian holonomies are a natural generaliza- 
tion of the Abelian Berry phases. We first assume that 
T is an iso- degenerate Hamiltonian family i.e., all the el- 
ements of T have the same degeneracy structure. This 
means that a generic Hamiltonian of T can be written as 
— J2h=i £ iW (A) where H;(A) denotes the pro- 



jector over the eigen-space W/(A) := span{|^ z Q (A))}™L 1 , 
with eigenvalues £;(A), whose dimension n\ is indepen- 
dent on the control parameter A. In order to preserve the 
R degeneracies ni we also assume that over M. there is 
no level-crossing i.e., I ^ V => £;(A) 7^ £j/(A), VA € M. In 
addition, we shall restrict to loops 7 in the control man- 
ifold i.e., maps 7: [0, T] i-> M such that 7(0) = j(T). 
These conditions in the dynamics of the system and in 
the control manipulations will facilitate the generation of 
holonomic unitaries. 

Let us state the main result on which the HQC 
relies. Consider a system with the above characteristics. 
When its control parameters are driven adiabatically i.e., 
slow with respect to any time-scale associated to the sys- 
tem dynamics, along a loop 7 in M. any initially prepared 
state \ipin) € H will be mapped after the period T onto 
the state 

IVw} = Ky \i> in ), U 7 = ffi^e** r Xl (7), (3) 

where, <pi :— dre;(A T ), is the dynamical phase 
whereas the matrices (7) 's represent the geometrical 
contributions. They are unitary mappings of TLi onto 
itself and they can be expressed by the following path 
ordered integrals 

r Al ( 7 ) :=Pexp I A t e U(m) , l=l,...,R. (4) 

J 'y 

These are the holonomies associated with the loop 7, 
and the adiabatic connection forms A\ . The latter have 
an explicit matrix form given by Ai = H;(A)<in;(A) = 
^^Ai^dXfj,, where Q analytically 

(^, M )^:=<C(A)|9/9A^|^(A)) (5) 

with (A^)^-! the local coordinates on M.. The connec- 
tion forms Aj's are nothing but the non-Abelian gauge 
potentials enabling the parallel transport |l5| over M. of 
vectors of the fiber 7i;. Result (||) is the non-Abelian 
generalization of the Berry phase connection presented 
first by Wilczek and Zee (1984) (see Appendix). Due to 
the decomposition of the evolution operator in (|^) into 
distinct evolutions for each eigen-space 7ii, we are able 
to restrict our study to a given degenerate eigen-space 
with fixed I. 

We shall present first an intuitive proof for deriving (^) 
and (|^) , aiming in clarifying the gauge structure interpre- 
tation of this adiabatic evolution and in providing a more 
physical insight. Without loss of generality we shall as- 
sume the family T to be iso-spectral. This implies that 
for any A € M it exists a unitary transformation IA{\) 
such that H{\) = U{\)H U(\)\ where H := H(X ). 
Upon dividing the time interval [0, T] into N equal seg- 
ments At, for Ui = U(j(X(ti))) one obtains the evolution 
operator in the form 
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U = Te~ 4 J* H ouH\)dt = T lim e -iJ2" =1 UiH ulAt 

N 

= T lim T[Uie~ iHoAt Ui (6) 

N—>oc A A 
i=l 

The third equality holds due to the smallness of the in- 
terval At in the limit of large N. The product LfjUi+i 
of two successive unitaries. gives rise to an infinitesi- 
mal rotation of the form Wlli+i » 1 + Ai ■ AAj, where 
(Ai)n = ll} £hf\ ■ The connection A has at time ti the 
components (Ai)^ with \i = 1, . . . , d. Hence the evolution 
operator @ becomes 

T lim U N ( t - iH N • At + V A\ ■ A\A U\. (7) 

For the case of a closed path the initial and the final 
transformations U\ and £//v are identical as they corre- 
spond to the same point of the control parameter mani- 
fold. With a reparametrization they may be taken to be 
equal to the identity transformation. Now we consider an 
initial state \tpi n ) belonging to an eigen-space Ho with as- 
sociated eigenvalue e.g., e$ — 0. Due to the time ordering 
symbol the actions on the state \tp(t)) of the Hamiltonian 
and of the connection A are alternated, hence in general 
we cannot separate them into two exponentials. On the 
other hand, if we demand adiabaticity, namely very slow 
exchange of energy during the process, this will keep the 
state within Ho, then at each time ti the state \ip(ti)) will 
remain in the £o = energy level. This allows to factor 
out in (0) the action of H, thus obtaining 




where A is projected into the subspace Ho- Notice that 
we replaced the time ordering with the path ordering P 
as the parameter of the integration at the last expression 
is the position on the loop 7. In this proof of the non- 
Abelian geometrical evolution it is clear how the holon- 
omy appears and which physical conditions enable its 
formation. In the same way we could have considered in 
addition to the equivalent transformations of the Hamil- 
tonian a multiplicative function £o(i) varying the energy 
eigenvalue. The results would be unaltered apart from 
the insertion of a dynamical phase. 

Let us now view some of the properties the holonomics 
have in terms of gauge reparametrization of the connec- 
tion and loop composition rules. In our context a lo- 
cal gauge transformation is the unitary transformation 
U(X) 1 > U(X)g(X), which does not change the Hamil- 
tonian operator Ho- Its action merely reparametrizes 
the variables of the control manifold. Taking into ac- 
count the properties gHo = Hog and = Tig we are 
able to obtain the transformation of the connection as 



A 1— ► Ag + g^ dg, (g S U{n)). It immediately follows 
that the holonomy transforms as Ta * g ^Ag- No- 
tice that in the new coordinates the state vectors \ip) 
i.e., the sections, become g^ This property makes it 
clear that the holonomy transformation has an intrinsic 
i.e., coordinate- free, meaning. Furthermore, the holon- 
omy has the following property in terms of the loops. 
We define (setting T = 1) the loop space at a given point 
Ao € M as 

L Xo := { 7 : [0, 1] 1 * M/j(0) = 7 (1) = A } 

over a point Aj £ M. Let us stress that, as far as the man- 
ifold M. is connected, the distinguished point Ao does not 
play any role. In this space we introduce a composition 
law for loops 

(72 • 7i) (*) = 6{\ ~ t) 7l (2i) + 9(t - i) 72 (2i - 1) (8) 

and a unity element 7 o(i) = Ao, t € [0, 1] moreover with 
7 _1 we shall denote the loop 1 1— * 7 (1 — t). 

The holonomy can be considered as a map Ta - L\ a t— > 
U(ni)j whose basic properties can be easily derived from 
eq. (ft: 

i) rU(72-7i) = r y i( 7 2) rU( 7 i); by composing loops in 
M. one obtains a unitary evolution that is the prod- 
uct of the evolutions associated with the individual 
loops, 

ii) T J 4(7o) = 1; staying at rest in the parameter space 
corresponds to no evolution at all, 

iii) T J 4(7 _1 ) = r^ 4 1 (7); in order to get the inverse 
holonomy one has to traverse the path 7 with re- 
versed orientation, 

iv) r ' a(ci 1 f) ~ -Tyi( 7 ), where tp is any diffeomorphism 
of [0, 1]; as long as adiabaticity holds the holonomy 
does not depend on the speed at which the path is 
traveled but just on the path geometry. 

From the properties listed above it is easy to show that 
the set Hol(A) := Ta(L\ ) is a subgroup of U(n). Such 
a subgroup is known as the holonomy group of the con- 
nection A. When the holonomy group coincides with the 
whole U (n) then the connection A is called irreducible. 
The notion of irreducibility plays a crucial role in HQC in 
that it corresponds to the computational notion of uni- 
versality [[Tj|. In order to evaluate if this condition is 
fulfilled by a given connection it is useful to consider the 
curvature 2-form F — dx^ /\dx v associated with 

the 1-form connection A whose components 

Fjxv =dnA v - d v A^ + [A„ , A v \ . (9) 

The relation of the curvature with irreducibility is given 
by the following statement | ^5[ : the linear span of the 
Fpv 's is the Lie algebra of the holonomy group. It follows 
in particular that when the F^s span the whole u(n) 
the connection is irreducible. 
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C. Holonomic Quantum Computation 

The unitary holonomies (Q) are the main ingredient 
of our approach to QC. From now on we shall con- 
sider a given subspace Hi (accordingly the label I will be 
dropped). Such a subspace, denoted by C, will represent 
our quantum code, whose elements will be the quantum 
information encoding codewords. The crucial remark |t]] 
is that when the connection is irreducible, for any cho- 
sen unitary transformation U over the code there exists 
a path 7 in M. such that ^^(7) — U\\ < e, with e arbi- 
trarily small. This means that any computation on the 
code C can be realized by adiabatically driving the control 
parameter configuration A along a suitable closed path 7. 

In particular we aim to constructing specific logical 
gates by moving along their corresponding loops. Ini- 
tially the degenerate states are prepared in to a "ground" 
state, interpreting the |0...0) state of m qubits. The 
statement of irreducibility of the connection A relates 
a particular unitary U with the loop ju over which the 
connection is integrated to give Ta^u) — U. Hence, 
there are loops in the control space such that the associ- 
ated holonomies give, for example, a one qubit Hadamard 
gate or a two qubit "controlled-not" gate. 

Let us emphasize the fact that one can perform uni- 
versal QC by only using quantum holonomies is remark- 
able. Indeed this kind of quantum evolutions is quite 
special, yet it contains in a sense the full computational 
power. On the other hand one has to pay the price given 
by the restriction of the computational space from H to 
its subspace C. Notice that, for the irreducibility prop- 
erty to hold, a necessary condition is clearly given by 
d(d — l)/2 > n 2 where d := dim.M. In particular this 
implies that for an exponentially large code C one has to 
be able to manipulate an exponentially large number of 
control parameters. 

Moreover like in any other scheme for QC, once the 
computation is completed a final state measurement is 
performed. To this aim it could be useful to lift the en- 
ergy degeneracy in order to be able address energetically 
the different codewords g] . This can be done by switch- 
ing on an external perturbation in a coherent fashion. 

We conclude this section by discussing the computa- 
tional complexity issue. The computational subspace C 
does not have in general a tensor product structure. This 
means that it cannot be viewed in a natural way as the 
state-space of a multi-partite system for which the notion 
of quantum entanglement makes sense. The latter, on the 
other hand, is known to be one of the crucial ingredients 
that provides to QC its additional power with respect to 
classical computation. It follows that, from this point of 
view, the scheme for HQC described so far is potentially 
incomplete. Indeed - as it will be illustrated later by ex- 
plicit examples - if N = dim C = 2 k i.e., we encode in 
C k qubits, then for obtaining with a multi-partite struc- 



ture a universal set of gates one needs 0{N) elementary 
holonomic loops. Thus in general one has an exponential 
slow-down in computational complexity. 

In Ref. H| we argued how one can in principle over- 
come such a drawback by focusing on a class of HQC 
models with a multi-partite structure given from the very 
beginning. The basic idea is simple: one considers an 
holonomic family T associated to a genuine multi-partite 
quantum system such that local (one- and two-qubit) 
gates can be performed by holonomies. Then from stan- 
dard universality results of QC p| stems that efficient 
quantum computations can be performed. An explicit 
example of the above strategy is formalized as follows 

Let us consider N qu-trits. The state space is then 
given by Hj ^ C 3 = span{|a)j / a = 0, 1, 2}. The holo- 
nomic (iso-spectral) family has the built-in local struc- 
ture T = {Hij(Xij)} where the local Hamiltonians ily 
have a non trivial actions only on the i-th and j-th factors 
of H. Moreover, Hij admits a four-dimensional degener- 
ate eigen-space dj := span{|a)i ® \f3)j / a,(3 — 0, 1} C 
Hi ®Hj = C 9 . If the Hij's allow for universal HQC over 
Cij then universal QC can be efficiently implemented over 
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C := Bpan{®£ 1 |a i ) j / a* = 0, 1} S (CT) 



III. THE CP N HOLONOMIC CONSTRUCTION 

In this section we shall consider a theoretical model 
where the holonomic ideas can be materialized. The CP™ 
model will be considered for which the Hamiltonian due 
to its degeneracy has such a symmetric structure as to 
allow the control manifold M. to be the n-dimensional 
complex projective space CP™. For quantum compu- 
tation we are interested in finding the particular loops 
which generate various holonomic gates and eventually 
constructing a complete set of universal gates. 



The path ordering prescription given in (A7) makes 
hard the explicit analytical evaluation of the holonomies. 
In order to tackle this problem we employ two pro- 
cedures. On the one hand we study loops restricted 
onto particular two dimensional sub-manifolds having 
easily calculated holonomies. This geometric restric- 
tion overcomes the difficulties connected with path or- 
dering, allowing the evaluation of a complete set of basic 
holonomies. On the other hand it is possible to com- 
pose a generic unitary operator with combinations of el- 
ements within this set. Eventually, by the loop compo- 
sition properties of the holonomies it is possible to find 
its corresponding composed loop. Even if the techniques 
presently known do not give the possibility to calculate 
the holonomy of the most general loop, we shall obtain 
families of loops and their corresponding holonomies from 
which any desirable group element may be constructed. 
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A. The CP n Model 

Consider the degenerate Hamiltonian Hq = e\n+l) (n+ 
1| acting on the state-space Ti. = C" +1 — span{| a) }",+:[. 
We shall take as the family T the whole orbit O(Hq) := 
{U HqW /IA £ U(n+ 1)} of Hq under the adjoint action 
of the unitary group U(n + 1). This orbit is isomorphic 
to the n-dimensional complex projective space 



O(H ) 



U(n + 1) ^SU(n + l) 



U(n) x U(l) 



U{n) 



CP' 



Each point, z, of the CP™ manifold corresponds to a 
unitary matrix U(z) — U\{zi)U2{z2)--- U n (z n ), where 
U a {z a ) = exp[G Q (z a )] with G a {z a ) = z a \a)(n + 1| - 
z a \n + and z a = # a e^=, for a = 1, ...,n. We shall 

assume in the following the set (6, <p) as real coordinates 
for CP™. The eigen-states of the rotated Hamiltonians 
are 

\a{0,4>)) :=U{6,4>)\a) = cos9 a \a) - 

n+l 

exp(—i(f) a )shi9 a ^ exp(i(j)j) sin6j jj^ cos# 7 |j) (10) 



j>a 



j>~/>a 



and 



|n + 1(0,0)) :=W(0,0)|n + l) = 

n+l 

exp(i(/)j) sinOj cos# 7 |j) (11) 



i=i 



7<i 



where we have defined n+ i := ir/2 and 4> n +i '■— 0. The 
first n ones © have zero eigenvalue while the last one 
( |TT| ) has eigenvalue e. Notice that for n = 1 the standard 
2-level model with the Abelian Berry phase is recovered. 



B. The Connection A and the Field Strength F 

By using definition (||) the components of the connec- 
tion A can be explicitly computed. Their particular form 
depends on the bundle of the degenerate spaces described 
in ( |lC)| ) and ([ll]). For CP™ A has 2n component as many 
as the dimensions of the manifold. These it(n)-valued 
connection components over CP™ are anti-hermitian ma- 
trices as dictated by (|J). In detail the only non-zero 
elements of the matrix A® 13 {(3 — 1,. .. ,n) are AJ^ for 



,{3- 1, given by 



A 13 - 



(a\U^U\(3) 



i(<t>s-<t>f>) 



sin V n, 



n 

P>1>a 



cosf 



(12) 



-A/ a which guarantees the anti- 



as well as A a ^ 

hcrmiticity. These are n components corresponding to 
the n ^-coordinates of CP™. The anti-hermitian matrix 
A<fip corresponding to the n ^-coordinates has non-zero 
elements for a — (3 and a > a given by 

i) 



-ij&o-M smQp S in( 



n • 

/3>7>a 



where we assumed Il/3> 7 >/3 

cos 7 — 1 , and for (3 > 

and a > a by 



Af a = {a\U^^-U\a 







' sin 8 a sin & sin 2 Op Y[ cos ®i Y\. cos 

/3>7>a /3>7>a 



Having the transformations U, which allow us to fix A, 
we are able to determine also the curvature and check 
the irrcducibility properties of the connection for the 
CP" model. By using the definition (||) and setting 
z a = + i z\ , one finds that at z = the components 
of the curvature are given by 



,(o) = n z [ 



dU a dU 



dz 



dz 



n 5 



z=0 



with a, (3 — 1, . . . , n and i, j = 0, 1. The relevant projec- 
tors are given by H z = W(z)IIW(z)', where n denotes the 
projectors over the first n degenerate eigenstates. Since 
dU a /dzi\ z=Q = i l (\a)(n+l\ + (-iy \n+l)(a\), one finds 



F 



— Al+3 



,j(0)=, 



{(-iy\a)(f3\-(-iy\P)(a\] . 



From this expression it follows that the components 
of F span the whole u(n) algebra. As remarked ear- 
lier this result does not depend on the specific point 
chosen, therefore the case considered is irreducible i.e., 
^a(L\ ) = U(n). Notice that in order to generate the 
loops in CP™ one needs to control 2n real parameters 
instead of the n 2 ones necessary for labeling a generic 
Hamiltonian. 

We can now open the way for applying these abstract 
constructions to quantum computing. In order to gener- 
ate a given quantum gate g € U(n) one has to determine 
a loop C g in M = CP" such that T A (C g ) = g. As the 
connection of the CP" model is irreducible, this is pos- 
sible for any group element g. Due to the non- Abelian 
character of the connection such an inverse problem is 
in general hard to solve. To tackle it we shall take ad- 
vantage of the composition properties described in the 
previous section. In particular one chooses specific fami- 
lies of loops {Ci}, that generate a complete set of easily 



G 



calculated holonomies, one can eventually construct any 
U (n) transformation. To this end first we consider the 
2-dimensional sub-manifolds in the 2n-dimensional space 
(0, cfi), spanned by two variables, (dp, (fig) or (9g, 9g), for 
specific values of j3 and (3. For these loops the line inte- 
gral, given by 



.4 



(A f>d6g + A x 0dXp) , 



(13) 



where A^ = 0^ or (fig, includes only two of the 2n com- 
ponents of the connection. Of course one cannot just 
simply calculate this line integral and then exponentiate 
it because of the path ordering procedure, which is nec- 
essary as the matrices A 6 f and A x e in general do not 
commute with each other. To overcome this difficulty 
we perform a second step for a further restriction of the 
loops {Ci}. From ([l2|) one checks that the parameters 
which define the position of the plane (9g, Xg), where the 
loop C lies, can be always chosen in such a way that the 
matrix A s f is identically zero. In particular, if one takes 
9i = 0, Vi 7^ (3, (3, the matrices A 6 ! 3 and A x v commute, 
so that one can calculate the integral and exponentiate 
it avoiding the path ordering problem. Of course the 






while for the (9p, (fip) components we have 

A 9 ? = 



and 



choice of the planes has to be such that the connection 
components lying on it do not give rise to a trivial holon- 
omy even if they commute with each other. This is indi- 
cated by the non-vanishing of the related field strength 
component Fg^x^. Another interpretation of the holon- 
omy, within this approach, is as the exponential of the 
flux of F, through the loop C. This definition is possi- 
ble once the problem has been "Abelianized" by having 
one of the two relevant components of the connection 
equal to zero. Application of the Stokes theorem then 
provides a natural way to evaluate the path integral, as 
§ c A x ed\p — f v , C s Fg^\pd9gd\p, where V(C) is the sur- 
face the loop C encloses on the (9g, A^-plane. 

In this framework, it is possible to identify four fami- 
lies of loops in such a way as to produce the basis of four 
matrices (the Pauli matrices and the identity) of all pos- 
sible two-by-two sub-matrices belonging to the algebra of 
U(2). With this approach one may restrict to a subspace 
of the degenerate space spanned by the states \/3) and 
|/3), ordered in such a way that (3 < (3. The relevant sets 
of coordinates are (9g, (fig) and (9a, (fig). Taking 9 L = 
for all i ^ (3, (3 one obtains the (9g, (fig) connection com- 
ponents 



—i sin 




% o 





sin^e^"^) 
'Sin^e- 1 ^-^ 



and 



A<t>B 



i sin 2 9g sin 2 tig 
-i sin 9g sin 9p cos 9pt 



-i{4> p -<t>p) 



-i sin 9 3 sin 9-g cos fl^e*^" - ^ 
—i sin 2 Op 



(14) 



(15) 



(16) 



With these four matrices we want to build a complete 
set of generators for the U(2) group. Specific choices 
of coordinate planes inside the four dimensional sub- 
manifold with coordinates {9g, (fig, 9p, tfip} shall provide 
those matrices. Note that this sub-manifold is locally 
isomorphic to CP 2 , and due to the irreducibility of its 



Ft 



sin 26^ 




relevant connection it is possible to produce the whole 
U(2) group. 

We can calculate the components of the field strength 
associated with the connection components ( 14|) , (fl5"|) , 
(filf). With a straightforward application of (|) we ob- 
tain the non-zero field strength components to be 



sin 29 r sin 9n 



|cos 9g sin 29 pe 



-i(<t>p-<t>p) 



— cos 9g sin 29 pe 



i{<t>l3-4>fj) 







Fe e = cos 9g 











= —i sin 9g cos 2 



e i(<t>p-<t>0) 
e -i(<f>i3-4>0) Q 



and 



F0.4,, = i 



sin 00 cos 2 9g sin 2 fl^e"^" 



F r h arhx = \ sin 9g cos 2 9g sin 20^ 



cos 2 9g sin 200 



e K4>0-<t>p) 
.g— Kip -dp) o 
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The field strength will be used in the following section 
to calculate the holonomies. 



C. The Holonomies 1^(7) 

Let us see how the restriction on various planes affects 
the relevant connection components (pl|),(|l5|),(16). It is 



logical to choose one of the plane coordinates to be the 
Op one as Ag^ = and hence it commutes with all others 
components, as required for implementing our strategy. 
The first choice is the plane (0p,(f>p), where the non- 
zero component of the connection is = — ism 2 Op. 
The second choice is the plane (Op,(j>p) for /3 > (3, with 
9g = 7r/2, giving a different connection with two non- 
zero elements, — ism 2 9p and Ag| = —i. Of course 
the latter element will give zero when integrated along a 
loop. For j3 > (3 both matrices are identically zero, and 
give rise to a trivial holonomy. With these two compo- 
nents and for appropriate loops one can obtain all possi- 
ble U(n) diagonal transformations. Indeed, for the loop 
C\ € (9p,4>p) we obtain 

r A (c 1 ) = eM-i\P)(P\x i ] 

Til denoting the area enclosed by C\ on a S 2 sphere 
with coordinates (20p,cf>p). This is exactly the Abelian 
Berry phase which could be produced if the state \0) were 
non-degenerate as it does not get mixed with the rest of 
the states. For Ci G (0p,4>p) we obtain analogously the 
holonomy 

r A (C 2 )=exp[i|/?)</?|E 2 ] . 

Recalling the constraint (3 < (3, we see that one can pro- 
duce n — 1 distinct holonomies from C2 type loops. 

To obtain the non-diagonal transformations one has 
to consider a loop on the (Op, Op) plane, with Oi = 
for all i 7^ (3, (3. Then, the only non- vanishing elements 



of the connection are A 



sin( 



= —At 



*PP " ' ",8/3 • 

By choosing further the (Op, Op) plane at the position 
<j>p = cf>p = the holonomy becomes, for the loop 
C3 G (Op^p)^^^ 



T A (C 3 )=eM-i(-iW)0\+A0)(0\)£ 3 ] 



while at cf>p = tt/2 and 
(O0,Op)< } , lj=7T /2^p=o we have 



0, i.e. 



r A (C 4 )=exp[-i(\f3)(P\ + \P)(p\)Z 4 ] , 

where E is the area on the sphere with coordinates 
(tt/2— Op, Op). The positive or negative sign in front of the 
area depends on the orientation of the surface enclosed 
by the loop C with respect to the orientation of the field 
strength F. Note that any loop C on the (0p,Xp) plane 



with the same enclosed area Ec (when mapped on the 
appropriate sphere) will give the same holonomy inde- 
pendent of its position and shape. These four holonomies 
restricted each time to a specific 2x2 sub-matrix generate 
all U(2) transformations. Hence, considering the inverse 
problem of obtaining a desired unitary from a holonomy 
we are able to choose from a whole family of loops of 
arbitrary shape and position. 

Finally, it is easy to check that in this way one can in- 
deed obtain U = exp[/i a T a ], where T a (a = 1, . . . , n 2 ) is 
a u(n) anti-hermitian generator and \i a an arbitrary real 
number. Therefore any element of U(n) can be obtained 
by controlling the 2n parameters labeling the points of 
CP". 

It is instructive to consider the form that the Hamil- 
tonian family T takes when restricted to the particular 
2-sub-manifolds. For the loop C\ (and similarly for C2) 
one finds 



Hi = -e/2B(20p,4>p)-a 

for B(0i,4>j) = ( sm Oi cos (j>j ;, sin Oi sin cf>j, cos0i) T , where 
the only non-zero elements are on the /3-th and (n + l)-th 
row and column. Hi generates an Abelian CP 1 phase in 
front of the state \f3) and the conjugate one in front of 
\n+ 1). On the other hand for the path C3 (and similarly 
for C4) we have 

H3 = eB(0p, Op)B(Op,0p) T , 

where the non-zero elements connect the states \f3), \p) 
and |n+ 1). In this Hamiltonian there is direct coupling 
between three states, giving rise to a non- Abelian inter- 
action. While H\ is easy to simulate in the laboratory 
with various experimental setups (spin 1/2 particle in a 
magnetic field, NMR, optical polarization, etc) it is yet 
an open challenge to construct the interaction dictated 
by the Hamiltonian H3. 



IV. APPLICATION TO THE OPTICAL 
HOLONOMIC SETUP 



In the previous sections we have developed the the- 
oretical background in order to make the non-Abelian 
geometrical phases useful into the quantum computing 
arena. Complex as it may be, such a construction of- 
fers various possibilities and advantages when applied to 
physical systems. In particular we shall resorts to quan- 
tum optics in order to make a physical application of 
HQC M, but also to clarify and resolve some theoretical 
issues discussed in the previous sections. To this aim we 
employ existing devices of quantum optics, such as dis- 
placing and squeezing devices and interferometers, acting 
on laser beams in a non-linear medium. 
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A. Displacers, Squeezers and Interferometers as 
holonomic devices 

In this subsection we shall consider the realistic imple- 
mentation of the holonomic computation in the frame of 
quantum optics. Even though the complete implemen- 
tation of the model presented here is likely to be exper- 
imentally a very challenging task, it is still remarkable 
that the necessary employed devices are realizable in the 
laboratory. Let us briefly present the setup. 

In order to perform holonomic computation with laser 
beams we shall consider the non-linear interaction Hamil- 
tonian produced by a Kerr medium 



Hi = hXn(n - 1) 
with n = a^a the number operator, 



a and a* being the 
usual bosonic annihilation and creation operators respec- 
tively, and X a constant proportional to the third or- 
der nonlinear susceptibility, > °f t ne medium. The 
degenerate eigenstates of Hi are the |0) and |1), where 
{If); u = 0,1, ...} denote the Fock basis of number eigen- 
states, n\v) — v\v). The degenerate space they span will 



• 1 • 



-iy 

cosh 2ri — e 1 sinh 2ri 



be the one qubit coding space. For the tensor product 
structure of I qubits we have to employ a set of I beams, 
providing the basis states \v\...v m ) = |z/i) ® ... ® \v m ) 
where v\ could be zero or one, for I = 1, m. The code 
can be written in this space of states. 

The iso-spectral transformations of the Hamiltonian, 
H(a) — U (a)HiW(cr), can be constructed by resorting to 
displacing and squeezing devices with unitaries -D(A) = 
exp(Aat — Aa) and S(/i) — exp(/.ta^ — p.a 2 ) respectively, 
as well as two mode displacing and squeezing devices: 
N(£) — exp(£a|a 2 — £aia 2 ) an d = cx P(C a i a 2 — 

£0102) respectively. In fact to obtain holonomies for one 
qubit gates we employ U(a) — D(X)S(fi) while for two 
qubits we take U(cr) — N(£)M((). In the first case the 
connection components are two-by-two matrices, where 
the null action is assumed on the rest of the tensor 
product sub-systems. For the coordinate decomposition 
A = x + iy and fi = r*i exp i9\ we have after some sim- 
ple algebra, similar to the one for CP 2 , the connection 
components 



-(cosh2ri — e lBl sinh2n) 
-iy 



ix z(cosh2ri + e 101 sinh2ri) 

i(cosh 2ri + e~ 101 sinh 2ri) ix 



. 1 








An 



1 

3 



■(cosh4ri — 1) 



The diagonal elements of the above matrices are also 
given easily by the equivalent Berry phases acquired by 
non-degenerate Fock states \v). For example the Abelian 
phases produced by a displacer are equal for all Fock 
states and given by 



(17) 



0Berry = f (V dx ~ ^dy) 



while the phases produced by a squeezer are given by 

<t>Berry = —J— j> (COsh^! - 1) d6 1 . 



Notice the complete similarity of the generators of the 
squeezing phases with Aq 1 as the non- Abelian matrix is 
diagonal, while in the displacing case only the diagonal 
elements are reproduced by (|l7|). 

On the other hand the parametric space of the two 
mode interferometers gives rise to the following connec- 
tion components 



A r2 — 



A r3 = 



as well as 



A 0n = 






































i6 2 
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-i0 2 
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3 



(2cosh 2 r 2 - 1) 



■ sinh 2r 2 



F (cosh2r 2 - 1) , 



and 



9 



-i9 3 






and 



i cosh 2r 2 sin 2r 3 




10 
0-10 




l sin r 3 



where £ = r2expi#2 and £ = r3expi#3. Here we define 
M := {(Ti\ th (Ti's being the real coordinates {x, y, r i: 9i} 
parameterizing all possible configurations with devices 
acting on the various laser beams. The first two com- 
ponents are enough for constructing holonomies repre- 
senting two qubit gates between any two qubits, which 
together with the one qubit rotations result into a uni- 
versal set of gates. It is straightforward to apply the 
conditions posed in the previous section for finding coor- 
dinate planes, with the relevant connection components 
commuting. As A ri is identically zero it can be com- 
bined with A x , A y and Ag 1 to produce U{2) holonomies, 
while A r2 and A r3 commute with each other leading to 
generate U{4) matrices. 

The holonomies obtained from the various loops are 
given in the following. The loop Cj € ( x > r i) g 1= o gi yes 

Fa (Cj) = exp-io-iE/ 

with £/ := f D{Ci) dxdr 1 2e~ 2ri . The loop C n € 
(2/>n) ei=0 gives 

T A (C n ) = cxp-i<r 2 X;// , 

with area £77 := § D( ^ C[ ^dydr\2e lri . The loop Cm S 
(riA) gives 

T A (C H i) = cxp -is 3 I]/// , 

with := $D(c in \d> r idQi sinh4ri. Cj, Cu and Cm 
produce U(2) rotations. In order to generate U(4) group 
elements we span the following configurations. The loop 



Civ S (r-2,r 3 ) 



,0 S lves 



iv) 



■ '-12-^ 

exp — ia 2 z->iv 



with Y*iv '■= J D ^ C[v ^dr2dr^2amh.2r2. The loop Cy G 



Tyi(Cy) = exp— ial 2 Y,\ 



with the area given by 
Above we have used 



S3 ■ = 



1 
3 



v dri dr-j, 2 sinh 2r 2 . 




-i 
i 
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while D(C P ) with p = I, V is the surface on the rele- 
vant sub-manifold (o{ , <jj ) of M. whose boundary is the 
path C p . The hyperbolic functions in these integrals stem 
out of the geometry of the su(l, 1) manifold associated 
with the relevant control sub-manifold. The T^C^'s thus 
generated belong either to the U(2) or C/(4) group and 
act on the one qubit space or on the space of the ten- 
sor product of two qubits, respectively. Considering the 
tensor product structure of our system these rotations 
represent in the 2 m space of m qubits respectively single 
qubit rotations and two qubit interactions, thus resulting 
into a universal set of logical gates. For example, Ta(Cv) 
with Sy = 7r/4 gives the following non-trivial two qubits 
gate 



V2 



V2 

1 -i 

-i 1 

V2 



while the holonomies produced by the loops Cj, Cu and 
Cm give a general one qubit rotation. Together they 
compose a universal set of transformations |2(J . 

The above description can be applied also to quantum 
systems with harmonic oscillator structure. For example, 
ion traps or atoms in a cavity are described by oscillat- 
ing modes (the ionic vibration or the cavity mode), which 
can be also displaced, squeezed and interfered with stan- 
dard experimental techniques [pT]-p4[. The degeneracy 
can be provided by the coupling of the oscillating modes 
with the internal energy levels of the atoms or by the 
structure of the product space of states of two modes 



B. SU(2) Interferometers and non-Abelian Stokes 
Theorem 



In what follows we discuss the possible use of SU(2) 
interferometer as a control devices pt| for producing 
holonomies. While relevant to quantum optical devices, 
such a scheme has an additional theoretical interest. 
Rather than distilling the Abelian sub-structure to avoid 
the path ordering problem we shall employ the non- 
Abelian Stokes theorem in order to rewrite the loop in- 
tegral of the connection as a surface integral of its field 
strength. The main advantage of this approach is the par- 
tial relaxation of the path ordering conditions, enabling 
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the holonomic calculation of non-commuting connection 
components. 

For ai and a 2 the annihilation operator of two different 
laser beams, consider the Hermitian operators 



Jx — \{a\a 2 + a\ai) , J y = -j{a\a 2 - a\a\) , 



and 



Jz 



N 



\{a[a\ - a\a 2 ) 



a\a\ + a\a 2 = n\ 



(18) 



n 2 



The operators ( |18| ) satisfy the commutation relations for 
the Lie algebra of SU(2); [J x ,J y ] — iJ z , [J y ,J z ] = i,J x , 
[Jz, Jx] — iJy The operator N, which is proportional to 
the free Hamiltonian of two laser beams, commutes with 
all of the J's. On the other hand, however, the Kerr 
Hamiltonian does not commute with the J's, allowing 
for the possibility that SU (2) interferometers be used as 
transformation controllers in view of the holonomic com- 
putation. 

From the above operators we obtain the unitaries, 
U x (a) = exp(iaJ x ), U v {(3) = exp(if3J y ) and U z (j) = 
exp(i-fj z ). For the degenerate state space of two laser 
beams spanned by \v1U2), we have from (^) and for 
U = U x (a)U y ((3)U z (^) the following connection compo- 
nents 





sin/3 cos/3e*'T 

cos/3e -17 — sin/3 







e 17 







and 
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-1 


















These components do not commute with each other, 
when projected on planes with non-trivial field strength. 
Hence, it is not possible to adopt again the method used 
previously to calculate the holonomies of paths in the 
three dimensional parameter space, (a,/3, 7). Instead, 
for this purpose we can employ the non-Abelian Stokes 
theorem p7] ]. The extra limitation, now, for the choice 
of the path comes from the constraint that, apart from 
being confined on a special two dimensional subspace, 
it has to have the shape of an orthogonal parallelogram 
with two sides lying along the coordinate axis. This will 



facilitate the extraction of an analytic result from the 
Stokes theorem. With the loop composition properties 
described earlier, it is possible to generalize these simple 
structures to any desirable loop. 

In order to state the non-Abelian Stokes theorem let 
us first present some preliminaries, where a few simpli- 
fications are introduced, as its general form will not be 
necessary in the present work. Consider the Wilson loop 
(holonomy), W — ~Pexp§ c A, of a rectangular loop C 
with the sides parallel to the coordinates (a, r) where 
a and r are a parameterization of the plane where the 
loop C lies. W is made out of the Wilson lines Wi 
for i = 1,...,4, as W = W4W3W2W1, where Wi corre- 
sponds to the i'th side of the rectangular ordered in anti- 
clockwise way. Define T _1 (er, r) = W4W3. Then, for F aT 
the field strength of the connection A on the plane (a, r) , 
W is given in terms of a surface integral 



W = P T e 



l (a,T)F^ T (a,T)T(a,T)d<TdT 



(19) 



where P T is the path ordering symbol with respect only 
to the r variable contrary to the usual path ordering sym- 
bol P which is defined with respect to both variables, a 
and r. Here F aT (a, r) = —d a A T + d r A a + [A a ,A T ]. 

The key point of this approach is that in (|l9|) the path 
ordering is only with respect to r so there is the possibil- 
ity to arrange the exponent in such a way that its inte- 
gration with respect to r is performed with a trivial path 
ordering procedure. Such possibility can be achieved in 
the present case by choosing the integration limits with 
respect to a so that they give commuting matrices as 
functions of r. In detail, we first choose a loop G\ on 
the plane (a, (3) positioned at 7 = 0. The relevant field 
strength component and the matrix T _1 (a, /?) are given 
by 



F a p — i(cos/3<7 3 — sin/^CTi) 



and 



T 1 = exp(i^a 2 ) exp — i — (sm(3&3 + cos/3cti) 
Hence we obtain 



T- x F af} T = 

— i(cos a sin 2(3 <j\ — cos a cos 2/3 03 + sin a a 2 ) (20) 

By taking a rectangle with the a-side equal to ir the only 
remaining term after the a integration is the one with a 2 , 
and hence the path ordering with respect to /3 does not 
pose any computational problem. As a result we have 
that for the closed rectangular loop C\ G (a, /3)-plane 
with coordinates {(0, 0), (a = 7r, 0), (a = n, (3), (0, /3)} we 
obtain the following unitary transformation 



T A (d) = expH2/3a^) 



(21) 
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With similar reasoning a rectangular loop C2 € 
(a, 7)-plane with coordinates {(0,0), (a = 7r,0),(a = 
7r,7),(0,7)} has the field strength component and the 
matrix T (a, 7) given by 



-i(cos7<T2 + sin 701 ) 



and 



T = exp (1 — 0-3 ) exp ~ (cos7<7i — sin 702) 
which give finally the holonomy 



T A (C 2 ) = exp(- l 2 7 a 3 12 ) . 

The matrix ct| 2 is defined similarly to a\* mm u 2 
Subsection IV A . Note that the coefficients in front of the 



matrices in the unitaries are areas on spheres spanned by 
the angles a and (3 or a and 7. This is consistent with 
the geometry of SU{2). 

Due to their structure the two matrices Ta(Ci) and 
^ a{C2) can produce any unitary transformation of one 
qubit encoded in the sub-space of states of two laser 
beams spanned by {| 01), 1 10)}. In other words, we need 
two laser beams to encode one qubit, contrary to previ- 
ous construction. The operation of two qubit interactions 
demands sophisticated optical devices as four mode in- 
terferometers, and we shall not discuss them here. 



C. Holonomies and Devices 

In the optical example studied above, a general state 
in the degenerate eigen-space of Hq = Hxerr is given 
as a linear combination of |0) and |1). The interaction 
Hamiltonian Hxerr describes fully our system when we 
move to the interaction picture. This transition affects 
only the parameters of the control devices which are re- 
defined. Under an iso-spectral, cyclic and adiabatic evo- 
lution of the Hamiltonian in the family J- ', the evolution 
operator acting on is given by the 2x2 sub-matrix 
in the upper left corner of the matrix 

U(0,T) = Texp-i / U{<r{t))H U\a{t))dt . (22) 
Jo 

This evolution takes place from time to time T and, to 
perform a closed path 7, we demand <r(0) = cr(T). The 
iso-spectral rotations of the Hamiltonian can be achieved 
with the use of the optical control devices described in 
the previous subsections. The degeneracy is provided 
by the propagation of the laser beams through the Kerr 
medium. Hence, when the laser beams are in the control 
devices, the degeneracy is not present. To restore degen- 
eracy along the adiabatic evolution one may employ the 
kicking method 0, [M now briefly recalled. 

Suppose that one is able to turn on and off a set 
of interactions (the "kicks") K := {U{a)} ae M on an 



ultra-fast time-scale with respect to the unperturbed part 
of the evolution. Let T — N At and to — 0, U + i = 
ti + At(i = 1, ... N — 1) be a partition of the time in- 
terval [0, T). The system evolution is as follows: at any 
time ti one kicks the system with the pulse U\ , j Ui where 
Ui := U{ai) is a unitary belonging to K (Uq = Un = !)■ 
If between the kicks the system evolution is unperturbed, 
one obtains the evolution along the loop 7 by 



N-l 



IL 



In the limit At i-> 0, N h-> 00, [N At = T] one gets U 1 -> 
Texp J T dtH(t) where H(t) := U{a{t)) H o U(a(0 . By 
making the function a(t) vary adiabatically, one can ap- 
proximate the holonomic evolution. The latter simulates 
the desired evolution ( f22|) by alternating the action of 
the degenerate Hamiltonian, Hq, and the control proce- 
dure with infinitesimal steps. We shall employ it here 
to demonstrate in a simple example the fidelity of the 
evolution operator approximated by the kicking method 
procedure and the relevant holonomy theoretically pre- 
dicted. 

Let us consider displacing devices D(X), performing a 
closed loop in their control parameters A. 




FIG. 1. The triangular (and polygonal) loop C on the 
complex plane of the displacing control parameters, A, ap- 
proximating the loop (circle). 

This is shown in Fig. |l|, where for simplicity the least 
possible number of displacing devices (three) for per- 
forming a closed loop has been drawn. When the state 
is at the edges of the polygon (triangle) no displace- 
ment takes place and the action of the Kerr medium 
is implied. Two displacing unitaries are combined as 
D(X)D(X') = exp (i3f(AA'))-D(A + A'). The physical pro- 
cess behind this is as follows. On the state lib) first 
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acts a displacing unitary D^(Xi), taking it to the point 
Ai. Then, the evolution operator of the Kerr Hamil- 
tonian, U(At) — exp(-iHoAt), acts for a time inter- 
val At = T/3. This effect is achieved by propagating 
the beam inside a Kerr medium. Then, the evolution 
Dt(A2)-D(Ai) is performed. This is obtained, with a sin- 
gle displacing device, given (up to an overall phase factor 
that will cancel at the end) by D(Xi — A2). After exiting 
the displacing device (we are at point A2) the beam en- 



ters a Kerr medium once more for time At and then the 
procedure is repeated until one comes back to the point 
Ai and the beam enters again the Kerr medium. Even- 
tually, the state is thus displaced by D(Ai). This loop 
may be transported to any other place of the control pa- 
rameter complex plane by acting at the beginning and at 
the end of this procedure with the appropriate displacing 
unitary (device). 



D 



D 



D 





X 




X 















FIG. 2. The setup with the alternating Kerr media and Displacing devices. 



The evolution operator is approximated by the follow- 
ing evolution operator 

£7(0, T) k, D{\i) (u(At; 0)U(At; A x + A 2 ) 
[/(A^AiMA^O)) £>t(Ai) , 



where [/(Ai; A) = D(X)U(At)D^ (A), A, = A m - A, and 
A4 = Ai. For integer N It is instructive to simulate nu- 
merically the above operator studying its dependence on 
the number of discrete steps used to approximate the 
continuous procedure. We start with five kicks. The in- 
volved parameters are taken to be T = 0.1 and hX = 1, 



with the radius of the circle equal to 1. The initial point 
is taken to be the origin of the complex plane rather than 
Ai, or in other words we do not perform the initial and 
final displacements by -D^(Ai) and D(Ai). In the table 
below are given the percentage deviations of the abso- 
lute values of the (0,0), (0,1), (1,0), (1,1) elements of the 
evolution operator {7(0, T) as functions of the number 
of displacers used to approximate the cyclic evolution. 
These are the relevant elements for the evolution of the 
states in the degenerate eigen-space describing a qubit. 
Numerical values were obtained for 5, 10, 20 and 26 dis- 
placers and the percentage is with respect to the value 
obtained with 100 displacers. 





5 


10 


20 


26 


00 


0.2419 % 


0.0595 % 


0.0149 % 


0.0099 % 


01 


0.9119 % 


0.2260 % 


0.0558 % 


0.0186 % 


10 


0.9119 % 


0.2260 % 


0.0558 % 


0.0186 % 


11 


1.6763 % 


0.4061 % 


0.0760 % 


0.0269 % 



We see that with 26 displacers the error is of the or- 
der of 1 in 10 4 acceptable for quantum computation with 
error correction. This provides an indication for the nec- 
essary number of devices needed in order to reproduce 
faithfully the holonomic adiabatic loop. 

In the polygon of Fig. ^ the displacement from point 
Ai to Ai+i is performed by the operating device with uni- 
tary action D(\i+i — A;). This device is naturally affected 
by an error in the position of the point Aj+i. As an over- 
all effect this will introduce an error in the area of the 
enclosed surface bounded by the loop C, as well as an 
error in the "matching" of the initial and final points Ai 
and Xn+i- The enclosed area is the parameter of the 
unitaries interpreted as logical gates, while the matching 
of the points Ai and Aat+i is required in order to apply 
the adiabatic theorem for the holonomic interpretation 
of the gates. In both cases the error is zero in the first 



order (i.e. our model is robust) with respect to the error 
introduced by spanning the loop C. 

Another appealing characteristic of the holonomic 
model is the following. Consider r^. to be the (kicking) 
time during which each displacer acts, T the overall pe- 
riod of the evolution and lo the characteristic frequency 
of our system with Hamiltonian Hq. Then, in order to 
apply the adiabatic theorem and the kick method the 
following inequality should to hold 

r fe < w -1 < T 

This condition is robust against small variations of the 
parameters Tk and T. 
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V. CONCLUSIONS 



VI. ACKNOWLEDGMENT 



In this paper we provided a detailed account of the 
so-called Holonomic approach to quantum computation 
(HQC). While in the standard "dynamical" picture the 
information processing is performed by a sequence (net- 
work) of logic gates obtained by turning on and off suit- 
able Hamiltonians in HQC the idea is to exploit the 
tools of non-abclian gauge theories to manipulate quan- 
tum information. The quantum codewords are realized 
by n-fold degenerate eigenstates of an iso-degenerate d- 
parametric family of Hamiltonians T = {H(X)}\. The 
quantum gates are then given by the holonomies U(-y) £ 
U(n) produced by moving along loops in the manifold M. 
of control parameters A. The holonomies represent the 
non-commutative generalization of the well-known Berry 
geometric phases and their existence is due to the non- 
trivial geometry of the set of the computational eigen- 
spaces thought of as a complex vector bundle over the 
control manifold M.. 

We focused on the case in which these loops are trav- 
eled in an adiabatic way. In this situation one can find 
explicit expression for the non-abelian (it(n)-valued) con- 
nection i.e., gauge potential, form A such that U(j) = 
P exp J A. The computational power of the given con- 
nection A is described in terms of the associated curva- 
ture form F : when d (d— l)/2 components of F span the 
whole u(n) then one can perform universal QC over the 
code. We thoroughly discussed several examples show- 
ing explicitly how to design loops, in the relevant control 
manifolds, in such a way as to get universal set of gates. 
We argued how, by using suitable Holonomic families T 
acting on multi-partite systems, one can achieve efficient 
computation taking advantage of quantum entanglement. 
We also discussed a potential experimental demonstra- 
tion of HQC in a quantum-optical set-up, using artificial 
Holonomic family T generated with an iterated "kicks" 
method. 

The HQC approach shows how many of the notions 
and techniques developed in (non-abelian) gauge theories 
can find a natural application and interpretation in the 
arena of quantum information processing. This rather 
unexpected connection between the tools used for the 
description of Nature at its fundamental level and the 
ideas of Information/Computation theory is conceptually 
quite intriguing. In a slogan HQC suggest that informa- 
tion, besides being "physical" it can also be "geometri- 
cal". On the more concrete side HQC, similarly to the 
other geometrical/topological approaches to QC recently 
emerged, shows inherent fault-tolerant issues. This lat- 
ter fact could make HQC appealing even from the imple- 
mentation point of view, in spite of the quite demanding 
requirements it involves. Individuation of potential can- 
didates for realization of HQC schemes in the lab is a 
major challenge for future research. 



The authors thank M. Rasetti and S. Chountasis for 
useful discussions and critical reading of the manuscript. 
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APPENDIX A: QUANTUM HOLONOMIES 

In this appendix we shall give a mathematical review 
of the basic formalism concerning quantum geometric 
phases and their non-abelian generalizations. Although 
the material discussed below is by now standard (see 
for example Ref. [ fl6[ on which our presentation strongly 
relies) it has been here introduced for the sake of self- 
completeness of the paper and for making it more acces- 
sible to readers from the field of Quantum Information. 
For the necessary (elementary) background of fiber bun- 
dles theory we refer the reader to the book [lq] . 



1. Abelian phases 

In quantum theory the physical states are represented 
by rays in a separable Hilbert space H. Mathematically 
this means that (pure) states are in a one-to-one cor- 
respondence with the elements of the projective space 
Pi/H)- The latter is defined as the quotient space of H 
with respect to the equivalence relation x ~ y <-* 3A G 
C — {0} I y — Xx, Vx, y EH. Alternatively one can con- 
sider the unit sphere in H, S°° := {\ip} / '\\\ip)\\ = 1}, 
factored by the U(l) action (IV'), e lB ) l_ * e t6 In this 
case one writes P(H) = S°°/U(l). 

The projection 



P(H) 



(Al) 



defines a U(l) principal bundle over the base space P(W) 
with total space 5°°. This fiber bundle is the natural set- 
ting in which abelian holonomies i.e., Berry phases, ap- 
pear. The fiber over the point |V , )(^I G P(?0 is given by 



Any other principal U(l) bundle can be written in 
terms of (Al) as follows. Suppose $ to be a (smooth) 
map from the parameter manifold M in P(TL), then 
one can construct the pull-back bundle §*S°°, with total 
space Uaeai ^*(A) an d projection 7T*: F^(\) —> A. A map 
A h-> \tpx) £ -f*(A) is a section of the bundle $*S°°. 

The Schrodinger equation id/dt\ip) — H deter- 
mines a temporal evolution 1 i— ► \ip(t)) in the total space 
5°° and via tt in P(H). Moving the other way round i.e., 
from an evolution in P(H) and obtaining an evolution in 
5°°, requires the introduction of a new key ingredient: a 
connection. 

The U(l) action over S°° defines the vertical direction 
at each point along the fiber. The connection w(l)-valucd 
1-form A allows to define the horizontal direction as well. 
Once this field of directions is given one can realize the 
horizontal lift of any curve in the base. 

To build the connection let us start by observing that 
S°° inherits from TL a natural hermitian structure. Such 
structure provides a natural mean for defining the hor- 
izontal directions at each point: the ones orthogonal to 
the fiber i.e., to \tp). Given the curve 1 i— > \ip(t)) <E S°° we 
decompose the tangent vector \ip) := d\ip(t))/dt as fol- 
lows \ip) = (ip\ip) \ip) + \h^). Where the horizontal com- 



ponent \hd)) satisfies the relation 



0. 



The connection can be evaluated explicitly by split- 
ting the operators d/dt according its vertical and hor- 
izontal components: d/dt = ad/89 + ^^B^ 
where ^ label the local chart coordinates (A M ) of the 
base manifold M.. The horizontal operators £> M are re- 
ferred to as the covariant derivatives, they are given 
by £> M := d/dX^ + A^d/89. The u(l)-valued 1-form 
A := A^ dX^ is the connection form. Apply- 
ing d/dt to \ip(t)) from the horizontal part one gets 



one obtains 



d/dX^ + A 
b\d/dX,\ 



M d/M\il>) = 



from which 
Since d/80\i>)=i \ip) 



(A2) 



A section s: M h-> 5°° /Ah of ( jAlj ) , under a local 
coordinate change, transforms as \ip)x l— * exp(i 8(X)) \ip)\ 
(gauge transformation). From this relation and eq. ( |A2| ) 
one can check that the connection form changes accord- 
ing to A^ i— > A^ — 88 /dXfj,. 

Let us consider now a loop 7: [0, T] 1— > P(Tt) in the 
base space. Using the connection form ( |A2| ) we can lift 
it to the total space S°°, let \ip(t)) be such an hori- 
zontal lift. One can write |V>(i)) = e 1 ^^ \tp(t)), where 
= (s o 7)(t), is a closed in S°° obtained com- 



posing 7 with a section. One has \ip(T)} 



\m), 



■/0e[0, 2tt)}. 



where (3 := f(T) — /(0). From the above relations and 
the horizontality condition (tp(t)\4>(t)) = 0, one gets (3 = 
£dtm)\iP(t)). Since \iP) = (9d/d9+j:Ld/dX^ ty), 
and 9(T) = 9(0) mod 2tt, by using Eq. (K2) one finds 
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p = E / dX ^ = / A - 



(A3) 



This expression of the [/(l)-holonomy is manifestly gauge 
invariant and it does not depend on the Hamiltonian: it 
is of a pure geometrical origin. 

In general the state vector is not horizontal {ip\ip) = 
—i (ip\ H (t )\ip) 7^ 0, in this case besides the geometrical 
term (A3) one has a dynamical contribution to the phase 

given by a := — J Q (^p(t)\H(t)\ip(t)) . At variance with 
the former this latter term depends on the Hamiltonian, 
moreover one can get rid of it by a gauge transformation 
\ip(t)} i-> U(t) \ip(t)), where 

U(t) :=exp[i f dr (i/>(t)\H(t)\1>(t))]. 
Jo 

By introducing the curvature 2-form F := dA = 
Yl^^d^ A d\ v , F^v = dAJ dX v - dA^/dX^, and 
using Stokes theorem, eq. ( [A3] ) can be rewritten as 
f3 := J E F where £ is any surface having 7([0, T} ) as 
a boundary. A non-vanishing F describes a non-trivial 
geometry of the bundle ( |Al[ ). One can generalize this 
construction to the more complex case in which the fiber 
is a n-dimensional complex space over which acts the 
group U(n). This non-Abelian situation will be now re- 
viewed. 



2. Non Abelian generalization 

Let us consider a family of Hamiltonians {H(X)}\<=m 
where M. is the control manifold. Wc shall consider loops 
7 in A4, and define H(t) := -ff 7 (t), t S [0, T]. In general 
H(X) — if 7 (t) has R different eigenvalues {ei}^ 1 with de- 
generacies {ni(X)}. We shall assume that no level cross- 
ing occur i.e., ni(X) — m. Let II; (A) denote the projec- 
tor over the eigen-space Tti(X) := span{| , 0"(A))}^ i =1 , of 
H(X). The spectral A-dependent resolution of the Hamil- 
tonians is then H(X) = J2i=i £ iW (A). The mapping 
A 1— > \ipf(X)} defines a section of the bundle 



U(N) 



U(N) 



U{N- ni ) U{ni)xU(N-m) 



(A4) 



where N := dimTt. The total (base) space of the U(ni)- 
principal bundle (A.4) is known as the Stiefel (Grass- 
mann) manifold and it is denoted by VN,m (Gn,ui ) ■ In 
the very same way as for the abelian case discussed above 
the hcrmitian structure over TL provides a natural notion 
of horizontality: tangent vectors are horizontal if they 



are orthogonal to the fiber. Notice that one recovers 
the abelian {TL = <C N ) case by setting m = 1, indeed 
G N ,i = CP^ 1 = P(H). 

The state vector evolves according the time-dependent 
Schrodinger equation idt\tp(t)) — H„m \ip(t)). In the case 
in which the loop 7's are traveled sufficiently slow one 
avoids transitions among different energy levels. In this 
adiabatic limit any initial preparation \ipo) G TL belong- 
ing to some energy eigen-space TLi will be mapped, after 
the period T, onto: \ip(T)) = U 1 |-0o) € TLi. For the sake 
of concreteness let us focus on the eigen-space associated 
with the eigenvalue si = 0. Also let {\4>a(t))} be the 
corresponding orthonormal basis at the instant t. 

Let n a (t) — J2f3U a p(t) ipp(t) be the solution of the 
Schrodinger equation with initial condition n a (0) = ?p a . 
By imposing, at each instant, the orthogonality i.e., hor- 
izontality, conditions (nplrja} = S a p one gets by differen- 
tiation 

S 

= u aS u; s + J2 u aS u; T (Vv \j> 6 ), (as) 

St 

from which follows (U' 1 U)p a = A a p where we have de- 
fined A a fj := (ipf3\4> a )- From the above relations it follows 
that 

U(t) = Texp [ drA(r) = P exp / A : (A6) 
Jo Jj 

where the ■u(n)-valued 1-form A = A^dX^ is given by 
{A^ap = (4> a \d / dX^ijjp) ■ Note that the matrix charac- 
ter of the connection A demands for the path ordering P 
to take place. For a reparametrization of the loop 7 in 
terms of a variable x £ [a,b] and for A(x) = A^-j^ it is 
formally defined by 



n=0 U - \ Ja ) 



with 



n 



b rti p%n-l 

dx\ I dx2--- / dxnA(xi)A(x2).--A(x n ) (A7) 

a •J a J a 



where Xk — a + ^(b — a). 
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